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Abstract 

We present a procedure of differential renormalization at the one loop level which 
avoids introducing unnecessary renormalization constants and automatically preserves 
abelian gauge invariance. The amplitudes are expressed in terms of a basis of singular 
functions. The local terms appearing in the renormalization of these functions are 
determined by requiring consistency with the propagator equation. Previous results in 
abelian theories, with and without super symmetry, are discussed in this context. 
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Differential regularization and renormalization (DR) Q was introduced as a renor- 
malization metliod in coordinate space compatible with gauge and chiral symmetry. In 
a series of papers this method has been further developed ^, ^ and successfully ap- 
plied to different theories [5-13]|^. However, it might be considered unsatisfactory the 
fact that Ward identities among renormalized Green functions are only satisfied when 
the different renormalization scales are conveniently adjusted. Instead, one would like 
that the gauge symmetry were automatically preserved, as occurs in dimensional reg- 
ularization and renormalization |T6| . 

In this letter we present a procedure to constrain the scales in DR at one loop 
while preserving abelian gauge invariance. This is done in two steps. First, each 
diagram is written in terms of a set of independent functions (with different number 
of propagators and/or different tensor structure). Second, the singular functions of 
this set are renormalized in a way which does not depend on the diagram where they 
appear. The local terms are fixed by the requirement that DR be compatible with 
the equation defining the propagator in the space of distributions. The propagator 
equation also allows to treat tadpole diagrams with the usual DR rules. In this manner 
one obtains renormalized Green functions which depend on just one arbitrary constant 
(the renormalization group scale) and, as we shall see, fulfil Ward identities in abelian 
gauge theories. The non-abelian case will be studied elsewhere. 

After describing the method, we discuss the renormalization of the one-loop vac- 
uum polarization in massive scalar QED, which is the simplest example requiring all 
the ingredients of the constrained procedure. The complete one-loop renormalization of 
this theory will be presented in Ref. ||T^. Then we review the one-loop Ward identities 
of massive QED and massless QED in an arbitrary gauge [^, the corresponding ABJ 
anomaly |8| and the evaluation of {g — 2)/ in supergravity, where supersymmetry is 
also preserved 

DR renormalizes diagrams by replacing singular expressions by derivatives of well- 
behaved distributions {differential reduction). These derivatives are understood in the 
sense of distribution theory, i.e., they are prescribed to act formally by parts on test 
functions {formal integration by parts). In practice, to carry out this programme one 
has to manipulate singular expressions. This gives rise to ambiguities which are usually 
taken care of by keeping arbitrary renormalization scales for different diagrams (or 
pieces of diagrams). The scales are adjusted at the end to enforce the Ward identities 
(which is equivalent to the addition of finite counterterms). In this letter we show that 
four rules are sufficient to formally manipulate and renormalize the singular expressions, 
^Different versions of difFerential renormalization can be found in nsll. 
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avoiding the introduction of unnecessary scales. The resulting renormalized amplitudes 
automatically satisfy the Ward identities. These rules are summarised as follows: 



1. Differential reduction, where we distinguish two cases: 

(a) Functions with singular behaviour worse than are reduced to derivatives 
of 'logarithmically' singular functions without introducing extra dimension- 
ful constants [^]. For example, 



(b) Logarithmically singular functions are written as derivatives of regular func- 
tions. At one loop we have the usual DR identity M 



1 1 logx^M^ 



which introduces a unique dimensionful constant (the renormalization group 
scale) . 

2. Formal integration by parts. In particular, 

[dF]^ = dF^ , (3) 

where F is an arbitrary function and R stands for renormalized. 

3. Delta function renormalization rule: 

[F{x, xi, ...,Xn)S{x - y)]^ = [F{x, xi, ...,Xn)]^6{x - y) . (4) 

4. The general validity of the propagator equation: 

F{x,Xi,...,Xn){0'' - m^)Amix) = F{x,xi,...,Xn){-5{x)) , (5) 



where Am{x) = i mKi{mx) ^ modified Bessel function |1£]. This is a 



valid mathematical identity between tempered distributions if F is well-behaved 
enough. This rule formally extends its range of applicability to an arbitrary 
function. 

The last rule will prove essential in our procedure. For instance, the 'engineering' 
tensor decomposition into trace and traceless parts is not compatible with it and will 
be modified by the addition of a finite local term. 
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To evaluate the diagrams we express them in terms of a set of basic functions 
using only algebraic manipulations (including four-dimensional Dirac algebra) and the 
Leibnitz rule for derivatives. To one loop these functions can be classified according 
to their number of propagators and their derivative structure. For the examples we 
discuss we need functions with one, two and three propagators: 

A = A{x)6{x) , (6) 
B[0] = A(x)O^A(x) , (7) 
T[0] = A{x)A{y)0-A{x-y) , (8) 

where A(x) = -^-^ is the massless propagator and O is a differential operator. In 
massive theories the following basic functions are also required: 

A = A{x)5{x) , (9) 
B[0] = A{x)0'-'A{x) , (10) 
T[0] = A(x)A(y)O^A(x-y) , (11) 

where A{x) = i^logx^m^ appears when the massive propagator is expanded in the 
mass m. Such expansion allows to properly separate pieces with different degree of 
singularity. The same type of functions also appear in massless theories if the photon 
propagator is written in a general Lorentz gauge. Note that A, A and B functions are 
singular, and B and T (T) are singular for n > 2 (n > 4), where n is the order of the 
differential operator O. 

We renormalize these basic functions using systematically rules 1 to 4. For exam- 
ple, 

T[a] = A{x)A{y)a^A{x-y) 
= -[Aix)f6{x-y) 

= -B[l]ix)5{x - y) ■ (12) 

T^p] = -B^[l](x)5{x-y) 

1 1 ^.logx^M^ , 

Our main observation is that the propagator equation (rule ^) can be further used to 
relate the different basic functions. Thus, by requiring that their renormalization be 
compatible with these relations, we shall completely fix the finite local terms (or scales) 
which appear in the differentially renormalized functions. Let us illustrate how to do 
this with one example, the renormalization of the basic function T[d^di,]. Using rule ^ 
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for a massless propagator, 

FnA{x) = -F5{x) , (14) 

we can write 

+ 2dyT[df,d^] + T[d^a] . (15) 

The last basic function on the right-hand side can be easily reduced to 



nd,o] = i,id;-dy)Tp] , (16) 



where rules ^ and |^ have been used. Now, we decompose the basic function T[df^du] 
into trace and traceless parts, adding an arbitrary (for the moment) local term to take 
into account the possible ambiguity introduced by this operation: 

T[d^d,] = ^ VT[n] + T[d^d, - i + ^b6ix)5{y)6^, . (17) 

The traceless part is finite because of the tensor structure and is not further renormal- 
ized. T[l] and T[5^] are also finite. On the other hand, the left-hand side of Eq. ( |l5|) 
is 

B[d^]{x)6{y) = A{x)d^^A{x)6{y) 



^cf^A{x)'5{y) 

1 

2 



-9^B[l](x)5(y) . (18) 



So, Eq. ( p^ ) reads for renormalized functions 

1 
2 



;9^B«[l](x)5(y) 



1 qxrri?r~-- „^.,„r„ „ 1 



+ I^bdyi6ix)6iy)) . (19) 

Since both members of this equation are finite, we can integrate on x using the inte- 
gration by parts prescription]^: 

= ldx{ayT[d,] + 2dyT[d,d^-^6,^a]) + ^^bdy6iy) 

^ A + lb)dy6{y) . (20) 



levr^M 2 



^Techniques for performing this sort of integrals can be found in Appendix A of Ref. ^ and Ap- 
pendix C of Ref. 
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Then this equation fixes 

(21) 

Note that the engineering trace-traceless decomposition, commonly used in the htera- 
ture of DR, is not compatible with the propagator equation in the case of logarithmic 
singularities. Hence, contraction of indexes does not commute with renormalization. 
Generically one must simplify all the tensor and Dirac structure before identifying the 
basic functions to be renormalized. The addition of the local term is equivalent to 
using a different mass scale M' in the renormalization of the T[n] coming from the 
trace-traceless decomposition, and then fixing log -^77 = b= — i. We prefer, however, 
to use the language of local terms to avoid confusion with the usual ad hoc adjustment 
of renormalization scales. 

With the same technique one can determine the renormalization of all the basic 
functions. In general, besides the massless propagator equation, Eq. (p^), one needs 

FnK{x) = FA{x) . (22) 

Both Eq. (|l^) and Eq. ( p2| ) are a consequence of the massive propagator equation, 
Eq. (^), and equivalent to the equation for the photon propagator in a general gauge, 

F(Vn - (1 - ^)d^,d,) A,p{x) = -F5^p6{x) , (23) 

where A^^(x) = 1/I67r2(5^^n + (a - 1)5^9;,) log x^^^^ 

The propagator equation can be further employed to 'separate' the tadpole func- 
tions into two-point functions, which can then be treated with the usual DR prescrip- 
tions: 

A = A{x)5{x) 
= -A{x)aA{x) 

= -Bp]. (24) 

In Table ^ we gather the renormalized expressions of the basic functions required 
in the applications below. In massive theories it is usually more convenient to work 
with compact expressions involving modified Bessel functions Q. The corresponding 
DR identities can be obtained by expanding the propagators in the mass parameter, 
using Table || and resumming the result. In practice, one uses recurrence relations 
among Bessel functions (see Appendix C of Ref. [|l2|) and then adds the necessary local 
terms to agree with Table ||. Table |2| collects the massive renormalization identities 
used in this paper. 
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= 

A« = -□B[l]+2a,B[5,] + i(2^ni2£|^ 

^ m - 4(4^^'-' — ^ — 

B^[n] = 

+ ii(4Sf(^m^--V°)'5W 
B^[n] = B^[i] 

T^P] = \j^^'^^S{x - y) 
Table 1: Renormalized expressions of basic functions. 
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B^[l] = ^{l{U-Am^ r''^^^-l^^^^-^ +^^\og^5{x)} 
+ 2vr2(logM^-l)5(x)] 

+ 27r2(log g + |)n<5(x) - 4vr2m2(l + Slog ^)5{x)]} 



Table 2: Renormalized expressions of massive basic functions, where = lS.m{x)5{x) 
and B„[0] = A™(x)0-A^(x). 

At this point any one-loop diagram can be renormalized: one just has to use the 
renormalized basic functions listed in the Tables. 

As a simple example which contains all the ingredients of the constrained pro- 
cedure let us consider in detail the vacuum polarization in massive scalar QED. The 
contributing diagrams are depicted in Fig. |l[ Using the Feynman rules given in Ref. ||l2| 
one gets 

n!i^(x) = -e2A„(x) S^Sz. A„(x) , (25) 
ng)(x) = -2eH^,yl^ra{x)5{x) , (26) 

which expressed in terms of basic functions read 



n|l](x) = -e^ABmid^d,] - d^dMi]} , (27) 
ng](x) = -2e25^^A^ . (28) 



Substituting the renormalized basic functions of Table ^, we obtain for each diagram 
4r(x) = -774w{(9.5.-V)[^C-4-^)(— "^"^"^'"^^"'"^ 



1 4 

+ m^Klimx) - Kf{mx))) + -7r2(log — - -)6{x)] 
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' 1 

W\A^; jAAAf WVAAAAAAAA; 

Diagram 1 Diagram 2 

Figure 1: One- loop diagrams contributing to the vacuum polarization of scalar QED. 

+ 5^,[27r^m\log^-mx)]} , (29) 
nJS''(^) = (i^2vrW(log^-l)V,5(x) . (30) 



The longitudinal terms cancel in the sum, yielding a transverse result: 



1 4 
+ m2(ir2(^:c) - Kfimx))) + -7r2(log — - -)5{x)] . (31) 

3 3 



In the same way, our method recovers previous DR results in abelian gauge theories, 
without the need to impose Ward identities a posteriori. Let us consider first the 
vacuum polarization in massive QED g. In ter ms of basic functions it reads 

U^,{x) = AemmH^, + ^Sf,,n - d^d,) B„[l] + 2B™[a^a,] - 5^,Bm[a]} (32) 



and the renormalized expression is 

R 



Ki^) = -v^id.d. - V)[^C - W)( -^o(mx)i^,(rn.) 



1 1 2 

- -m\Ki{mx) - Kl{mx))) + -vr2(log^ + -)5{x)\ , (33) 

which is again transverse. 

In supersymmetric QED the vacuum polarization is the sum of the spinor QED 
diagram and twice (two complex scalars for each Dirac spinor) the scalar QED diagrams. 
In terms of basic functions this gives directly a transverse result depending on one basic 
function only: 

n^,(x) = -2e^{d^d, - V°)Bn^[l] . (34) 
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In this case one could use an engineering trace-traceless decomposition in renormalizing 
each diagram. The gauge-non-invariant terms vanish in the total sum due to supersym- 
metry cancellations. The complete result would thus be the same as the one obtained 
from the renormalized functions of Table 0: 



M2 



+ 2tt^ log ^5{x)] . (35) 

Next we consider the QED vertex Ward identity between the electron self-energy 
and the electron-electron-photon vertex in an arbitrary Lorentz gauge, which was stud- 
ied in Ref. In this case and the next one the masses play no relevant role, as far as 
renormalization is concerned, so we consider massless electrons for simplicity. In order 
to respect the Ward identity for both the a-dependent and a-independent pieces (where 
a is the gauge parameter), the authors of Ref. |§ had to impose two relations among 
different scales: 

log^ = i> (36) 



Mv 4 
Me 
1^ 



A^logff = 3, (37) 



where My and Ms and appear in the vertex and in the two pieces of the electron 
self-energy, respectively. With the constrained method we find 

^ 9rl 1 o> log X^M^ 

+ i»-^ii\-(iiTf2»°^-^^^ + -[t^»m (38) 



for the electron self-energy, and 



1 



V^^{x, y) = ie3{-27,7^7,(9:9,^T[n] + S^f [5,a] - 5,^f [5„a]) + 4jaT[dad, - ^6a^O] 

1 1 logX^M^ ,11 e/ N 

+ (a - I)bp7a7^.7bla{d'ad!f[dpd^] + d',f[dbdpd^] - dlfidadpd^]) 

1 1 logx^M^^, , , 

(39) 

for the vertex. These renormalized amplitudes automatically satisfy the Ward identity 
{d; + dy)V^^ix,y) = ieJ:''{x-ymx)-6{y)) , (40) 
as can be seen by integrating on y. 



in 



The chiral triangle anomaly in QED was discussed in Refs. [Q, The renormal- 
ized triangle diagram depended on the relation between two renormalization scales. 
By adequately choosing these scales one could respect either vector current or axial 
current conservation, but not both. Imposing conservation of the vector current, the 
correct value of the axial anomaly resulted. In contrast, in the constrained DR method 
everything is determined and a non-ambiguous result is obtained: 

+ ^^,uxa{d:-dymx)5{y))] , (41) 

where the index A corresponds to the axial vertex. Then (see appendix B of Ref. |]l|) 

dlTl^.x{x.y) = 0, (42) 
dlTf^ux{x,y) = 0, (43) 

o3 



+ 9f)r^^,,(x,y) = '-^e,,x,dldy{5{x)5{y)) , (44) 

so the vector Ward identities are directly preserved while the axial one is broken, giving 
the known result for the anomaly. 

Finally, let us comment briefly on the calculation of the {g — 2)i in unbroken 



supergravity performed in Ref. [12|. There, the symmetry to be preserved was super- 



symmetry, which implies a vanishing anomalous magnetic moment |20|. This result 
was obtained thanks to the use of the propagator equation to explicitly relate dia- 
grams with different topology. Then, only one type of singular basic function, T[n], 
appeared. Although engineering trace-traceless decompositions were performed at in- 
termediate steps, this did not affect the total sum because the extra local terms cancel, 
as occurs in the vacuum polarization in supersymmetric QED discussed above. Using 
the renormalized basic functions in Tables [l| and |2|, {g — 2)i also vanishes, although the 
contribution of each diagram is different, as is the total graviton contribution. 

Summarizing, we have proposed a procedure of differential renormalization to 
one loop which only introduces a single renormalization scale. We have verified that 
the renormalized amplitudes so obtained automatically satisfy the Ward identities 
of abelian gauge symmetry in known examples, that the chiral anomaly is correctly 
treated, and that supersymmetry is preserved in a relatively complex calculation. In 
practice, one just needs to use the renormalized functions of Tables |l] and |^. 

In principle, the method could be generalized to higher loops. However this is not 
straightforward. New more complicated functions emerge, which could be tackled with 
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the systematic method of Ref. [Q. Still, one should properly constrain the local terms 
using rules 1 to 4 or some consistent extension of them. 

We thank J. Collins and R. Stora for discussions. This work has been supported 
by CICYT, contract number AEN96-1672, and by Junta de Andalucia, FQMIOI. RMT 
and MPV thank MEC for financial support. 
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